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Abstract. Coherence siniphccs are generic topological correlation- function defects 
supported by a hierarchy of coherence functions. We classify coherence simplices 
based on their topology and discuss their structure and dynamics, together with their 
relevance to several physical systems. 
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1. Introduction 

The quantized vortex is an archetypal topological defect which arises in a variety of 
physical systems including superfiuids, superconductors, optical speckle fields, coherent 
fields in the focal volumes of lenses, coherent optical fields scattered from sharp edges, 
eigenmodes of waveguides and angular momentum eigenstates of the hydrogenic atom 
[H [21 [3] . The structure and dynamics of these conventional quantized vortices bear a 
close similarity to classical eddies, while enabling some unique material properties such 
as crystallization of magnetic flux in type-II superconductors and quantized circulation 
in rotating superfiuids 

Let us focus on this quantization of circulation, restricting the discussion to complex 
scalar fields such as the wavefunction of Schrodinger wave mechanics or the complex 
analytic signal associated with a scalar monochromatic electromagnetic wave. Quantized 
vortices in such fields are very generic, as was emphasized by Dirac in a visionary paper 
from the 1930s [5] . When studying vortices from a topological perspective one can ignore 
the details of the particular differential equation obeyed by a given complex-valued field 
ip, including whether the said equation is linear or non-linear, instead relying on the 
very general physical property that ip must be a continuous and single- valued function of 
position and time. Using this assumption alone, Dirac was able to show that quantized 
vortices will in general exist in complex fields, and that they are stable with respect to 
perturbation. The essential idea, here, is that while the complex field is single valued, 
its phase need not necessarily be single valued; hence the integral of the phase around 
a closed loop need not vanish, and may in fact be an integer multiple m of 2tt radians. 
While we explore this point in further detail in Section 2, the key point to note here is 
that non-zero m heralds the presence of a quantized vortex. 

As the previous list of examples shows, quantized vortices in complex fields have 
long been present in physics. Examples include Schrodinger's angular-momentum 
eigenstates of the hydrogenic atom, Wolther's 1950 treatment of vortices in the context 
of the Goos-Hanchen effect, and the 1952 study by Braunbek and Laukien investigating 
the exact solution due to Sommerfeld, for the diffraction of a plane electromagnetic 
wave from an infinite half plane [B]. Nye and Berry's work on dislocations in wave 
trains [7] was pivotal in bringing the study of quantized vortices to the attention of the 
scalar optics community. The study of such optical vortices has now blossomed into 
the area of singular optics, with several key reviews including Berry P|, Nye [5], Soskin 
and Vasnetsov [10], and Dennis et al. |6]. The nodal lines which thread vortex cores 
can form closed loops or extend to infinity [S]; they can terminate at points where the 
potential is discontinuous; and they can even form knots [HI [121 1131 [IH [15] . 

Now, partially coherent optical fields or mixed-state wave functions are often 
studied via correlation functions such as the mutual intensity or the cross-spectral 
density (for partially coherent optical fields) [Ml [17] and the density matrix (for mixed- 
state wavefunctions) [TS]. These correlation functions are very often continuous single- 
valued complex functions of pairs of spatial coordinates. If one recalls Dirac's argument 
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for the existence of vortices in complex functions, one has the logical possibility of 
vortex-type structures in field correlation functions. Such vortical structures, in two- 
point correlation functions of both classical and quantum-mechanical complex fields, are 
known as coherence vortices. 

Accordingly, quantized vorticity has emerged as a topic of interest in the context of 



optical and matter field coherence [iniEOlEIlESlEHlEaESlESlEZlEHl^ For a 



recent review of coherence vortices, see Gbur and Visser [52] ■ Here the term 'coherence 
vortex' has been coined to describe phase singularities or nodal lines in the cross spectral 
density matrix, and related coherence functions such as the spectral degree of coherence 
[2U] . Coherence vortices may undergo topological reactions such as sphtting, fusing and 
pairwise creation/annihilation [28]. They may be spontaneously formed in a Young-type 
three-pinhole interferometer [27], and via Mie scattering from one or several spheres 
[29l 130] . Coherence vortices have also been studied in the context of hnear optical 
imaging systems and the focal volume of aberrated lenses [26]. Recently, it was shown 
that coherence vortices emerge even in systems with only one spatial dimension where 
conventional vortices are manifestly absent [31] . Coherence vortices have been observed 
experimentally in optical fields [211 [231 121] ■ 

All of the previously cited work refers to screw- type and/or edge- type topological 
defects in two-point coherence functions. The desire to investigate both higher-order 
correlation functions and more complex topological correlation defects has motivated 
the present investigation. Here we develop the notion of a 'coherence simplex' as a 
generalization of the concept of a coherence vortex, thereby categorizing topological 
defects emerging in generic p-point complex correlation functions. 

We close this introduction with an overview of the remainder of the paper. Section 
2 gives a brief background on conventional quantized vortices as screw-type phase 
singularities in continuous single-valued complex functions. In Sec. 3 we define a 
coherence simplex as a generalized form of coherence vortex, and in Sec. 4 we categorize 
a hierarchy of coherence simplices. Finally, in Sec. 5 we discuss the structure and 
dynamics of coherence simplices, followed by a discussion and concluding remarks in 
Sec. 6. 

2. Conventional quantized vortex 

To be self-contained and to facilitate further discussion, we briefly review some 
mathematical properties of conventional quantized vortices. A generic continuous 
single- valued complex scalar function / of two real variables (x, y) can, without loss 
of generality, be cast in the Madelung form 



where 0{x,y) is a real valued phase function [33] • If in the vicinity of some point (xq, yo) 
the phase function has the form 



f{x,y) = \f{x,y)\e'''^^'y\ 



(1) 
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where m is an integer and B{x,y) is any real function that is analytic in the vicinity 
of {xo,yo), then the phase singularity at (x = Xo,y = yo) is called a point vortex. 
Physically, we can view B{x,y) as a continuous background deformation, upon which 
sits the screw-type phase dislocation given by term proportional to m. The amplitude 
|/(xo, yo) I must vanish at the location of the phase singularity to ensure single valuedness 
of the complex function. In three-dimensional systems such zero-amplitude phase 
singularities are nodal lines or line vortices which may close on themselves forming 
loops or vortex rings P El |9l |3ll [35l [36] or even vortex knots [HI |12l |13l HU ITS] . 

In superfluids, the gradient of the phase function V arg['?/'(r)] of the complex order 
parameter field ip{r) describing the system can be associated with the velocity field fs(r) 
of the superfluid. Since the curl of the gradient V x Vv4(r) of any vector field ^(r) is 
identically zero unless the field contains singularities, applying Stokes' law directly leads 
to the quantization of circulation of the superflow 



where m is an integer and k is the quantum of circulation. In helium II and atomic 
Bose-Einstein condensates it is the circulation of atoms that is quantized [U ESI 137] . 
In superconductors magnetic flux is quantized and vortices form due to the motion of 
Cooper pairs [21 |38l [39], and in coherent optical fields the photons can form optical 
vortices [HI |36] . It is also possible to convert optical to matterwave vortices and vice 
versa 

From a mathematical perspective, any single-valued continuous complex function 
of at least two real variables may possess screw-type phase singularities of this kind [5] 
and hence vortices should be expected whenever a model of a physical system involves 
continuous complex functions. It is particularly interesting to investigate the vortex 
properties of complex p-point correlation functions of quantum fields, to which topic we 
now turn. 

3. Coherence simplices 

The previous section has restricted consideration to vortices in field functions, considered 
as a function of position. Other control parameters such as time may be present, 
but this does not change the fact that the vortices of Sec. 2 'swirl' in physical 
space; in order to swirl in physical space, either two dimensions {x,y), or three 
dimensions {x,y,z), are required. Yet, as was emphasized in the first two sections 
of this paper, the existence of vortices follows directly from one having a continuous 
single-valued complex function. This raises the logical possibility that complex two- 
point correlation functions, which are continuous functions of two spatial coordinates, 
{xi,X2), {xi,yi,X2,y2) or {xi,yi, zi,X2,y2, Z2), together with any other relevant control 
parameters such as time coordinates, may also possess vortical structures in the phase 
of the said correlation function. However, these vortical correlation-function phases — 



termed 'coherence vortices' [IHl[2ni[2Il[221[231l23[2Sl[2Sl[2Zl[2Hl 132]— will 'swirl' in 




(3) 
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spaces of high-dimensionality, rather than in physical space. Moreover, these topological 
correlation-function phase maps are not restricted to two-point correlation functions, but 
may evidently be generalized to p-point correlation functions. Such structures, termed 
coherence simplices, are the topic of this section. 

Let us adopt the language of second quantization and introduce Heisenberg 
field operators ip'^{r,t), '?/'(r, t) which respectively create and annihilate an excitation of 
the field, such as a particle, at space-time point (r, t). Furthermore, a bosonic field 
ip{r,t) obeys the canonical commutation relations 

[^P{r,t)J\r',t')]=6{r-r')6{t-t'), (4) 

with similar relations holding for fermionic fields with the commutators replaced by 
anticommutators [13] . 

Consider the p-th order field correlation [HI |371 HJl |13] (cf. many-body 
density matrix, Glauber quantum coherence functions. Green functions or Feynman 
propagators) 

g^P\xi,X2, . ..X2p) = (^"^(ri,ti)V^1'(r2,t2) • • • (5) 

where Xj refers to both space and time coordinates and the angular brackets denote 
quantum statistical average. We augment this by further defining 

(V^(r))=^?W(xo) (6) 
which may acquire finite values through the spontaneous symmetry breaking mechanism 

With these prerequisites, we may define a p-th order vortical coherence simplex as 
a generalized vortex in the complex valued correlation function g^'P\xi^ X2, ■ ■ ■ X2p)- For 
such vortical coherence simplices, coherence circulation is quantized via 

y V aTg[g^P\xi, X2, . . . X2p)] ■ dhpo = rriKp, (7) 

where dl2pD is a line element along an oriented closed curve O in a 2]9Z}-dimensional 
space, where D is the physical dimension of the space, and the constant Kp is a quantum 
of coherence circulation. 

As in the case of conventional vortices [3511111 HS], coherence simplices with winding 
number |m| > 1 are likely to be topologically unstable with respect to perturbations 
(cf. and split into a number of lower winding-number excitations. By the term 

coherence simplex we generically refer to any kind of topological defect in any correlation 
function and thus include e.g. soliton-like structures. 

In particular, the p = VCS reduces to a conventional quantized vortex in two- 
or three-dimensional space. It is immediately clear that if the field ip{r) contains a 
conventional vortex, it may also show up as a coherence simplex in the associated 
correlation functions (see [20l EH [23l El])- However, coherence simplices and VCSs in 
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Figure 1. Hierarchies of vortical coherence simplices in the space spanned by spatial 
dimension versus the order p of the coherence function. The lower left corner of 
each matrix element shows the codimension Cy of the vortical coherence simplex and 
the diagrams depict the simplex structure of the underlying coherence function. The 
numbers inside the bullets denote the physical dimension of the space. 



particular, may also appear when the underlying field itself possesses no conventional 
vortices. In fact, VCSs emerge even when ip{r) has only one spatial dimension and 
therefore could not exhibit conventional vorticity even in principle |31j . 

In the generalized core, or hypercore, of the VCS g^^'' vanishes and the phase 
winds an integer m times 2tt around it. While a conventional material vortex is a 
zero probability of finding a particle in some point which defines the location of the 
vortex core, the hypercore of a VCS is a zero probability of finding coherence between 
two (or more) space-time points. The coherence simplex can also be described as a 
topologically unavoidable loss of coherence between sets of space-time coordinates due 
to quantization of circulation of the fiow of coherence. In particular, consider a line 
integral, ([T]), which yields a nonzero coherence circulation winding number m. Even if 
the quantum field maintains full coherence at every point along and outside the closed 
curve O, the topology of the problem dictates there to be at least one generalized point 
inside O at which the system is fully incoherent, as is further clarified in Sec. 5. 

4. Hierarchies of coherence simplices 

Figure [T] illustrates the simplest topologies in the infinite hierarchy of possible coherence 
simplices. The dimensionality of the physical space is depicted on the vertical axis and 
the horizontal axis labels the order p of the correlation function g^'''\xi, X2 ■ ■ ■X2p) in 
which the coherence simplex is embedded (see e.g. ([6])). The dimension M of the 
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embedding space of a VCS is M = D for p = and M = 2p x D otherwise, where the 
dimension of the physical space is denoted hj D. A generic coherence simplex has a 
codimension 

C = M -N (8) 

where is the number of degrees of freedom consumed by the topological defect 
type. Vortical coherence simplices categorized in figure [T] therefore have a codimension 
= M — 2 which is denoted in the lower left corner of each pane. The diagram in 
each pane illustrates the simplex structure of the correlation function underlying the 
topological defect. Solitonic coherence simplices have a codimension Cg = M — 1 and, 
unlike coherence vortices [IHlEn], exist also in the D = l,p = case. However, solitons 
are not in general topologically protected in the way vortices are and therefore their 
stabihty with respect to perturbations |9] is determined by energetic considerations 
rather than their topology. Higher order N > 2 coherence simplices can be categorized 
using similar diagrams to that shown in figure [Tj 

The first (shaded) column in figure [l] corresponds to the conventional vortices, 
which are manifestly absent in systems with one spatial dimension. In two-dimensional 
systems these conventional vortices are point charges whereas in the three-dimensional 
case they correspond to line-like vortex filaments. The second column in figure [T] contains 
all possible VCSs supported by the two-point correlation functions. In particular, in ID 
these correspond to coherence vortices [5T] and in 2D they are sheet-like structures. 

Coherence simplices can be constrained to lower dimensional spaces revealing 
a correspondingly simpler topological structure. For example, slicing through a 
conventional 3D vortex line results in a 2D point vortex. Similarly, fixing three of 
the six coordinates {x,y, z,x' ,y' , z') of a first order p = 1 VCS in a physical space 
with three spatial dimensions D = 3 reveals the VCSs as a one-dimensional nodal line 
object, whereas if all of the coordinates are left free the VCS is a four dimensional object. 
Most of the coherence simplices discussed in the literature so far [TSl 1201 EH 123 12S1 ED] 
correspond to such reduced-dimensional or constrained VCSs. Curiously, as is evident 
in figure [TJ a line-like Cy = 1 VCS only exists for p > in the presence of constraints 
fixing one or more degrees of freedom but never as a free object. On the other hand, 
under suitable constraints any VCS with finite codimension Cy > can be constrained 
to appear as a line-like coherence vortex. 

5. Structure and dynamics of coherence simplices 

The zeroth-order simple vortical coherence simplices are simply conventional vortices 
with vanishing field amplitude in the vicinity of a screw-type phase singularity. In other 
words they correspond to a hole around which the entities described by the field, such 
as atoms or photons, circulate. In contrast, the higher order {p > 0) VCSs typically 
have finite probability of finding corpuscles inside their hypercores, despite the fact that 
the corresponding coherence function vanishes there. As in the case of conventional 
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Figure 2. Vector coherence simplices in a two-member (a,b) and three-member (c,d) 
ensembles of two-state wavefunctions defined in the text. The four different spin- 
spin correlations are denoted in the upper corners of each frame in each subfigure. 
Subfigures (a) and (c) show the coherence density and (b) and (d) show the phase of 
the coherence function. Two different kinds of vector coherence defects are circled in 
frames (b) and (c). The constant a is an arbitrary spatial length scale. 



vortices, due to energetic considerations multi-quantum {\m\ > 1) VCSs are likely to 
split into multiple lower order VCSs and multiple VCSs in equilibrium may arrange in 
a VCS lattice whose structure may be that of an Abrikosov-type triangle shape, square, 
etc. depending on the details of the defect topology and interactions. 

If the complex field supporting the VCSs has multiple internal degrees 
of freedom, the field operator becomes a multicomponent spinor '?/'J.(r,t) = 
('?/'|(r, t), ?/'|(r, t), . . . '?/'J.(r, t))-'", where T denotes transpose and a is a spin index, and 
hence the coherence functions have a tensor structure gl^},{xi,X2 ■ ■ ■ X2p). Through 
this construction such coherence functions will in general admit coherence singularities 
analogous to monopoles, skyrmions, fractional vortices, hedgehogs, sheets, textures, 
knots etc. (16| HTJ HHl IH]. In such spinor fields with vector order parameter, continuous 
coherence simplices for which even the coherence density does not have to vanish in the 
spinor VCS hypercores are possible. 

Figure 2 is an illustration of vector coherence defects produced by the mixture of 
(a,b) two and (c,d) three traveling two-component Gaussian wave packets of the form 



) / Hn, where is a normalization constant. To 
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produce figure 2 (a,b) we used dimensionless values xi = —1, fc| = fcj = 0.5, yUi = 1 and 
X2 = l,kl = = -0A,fi2 = 1 and for (c,d) Xi = 3.0, kl = -2.0, A;| = 2.0, /^i = 3.0, 
X2 = 0, kl = k^ = —1.0, /i2 = 2.0 and x^ = —3.0, ^3 = 1.5, ^3 = —1.5, fi2 = 1-5. 

Frames (a) and (c) sliow tlie resulting coherence densities \9^^}>{x,x')\ and (b) and 
(d) are the corresponding phase functions. The circles in (b) and (c), respectively mark 
the location of a winding number w = {1, 1, 1, 1} and w = {0, —1, 1, 0} defects where the 
generalized winding number w = {rrin, mi2, . . . mg-o-'} is expressed as vector of winding 
numbers corresponding to the different spin-spin correlation functions. 

Coherence simplices are generically dynamical objects flowing with the coherence 
current [23]. When perturbed, equilibrium coherence simplices are expected to 
reveal elementary excitation modes such as helical Kelvin waves [371 EHl EI] and 
Tkachenko shear waves [H21 E31 [Ml ES] and their VCS displacement- wave generalizations. 
In this context, figure 3 of Marasinghe et al. [30| displays a clear signature of 
oscillatory dynamics, similar to the Crow instability of antiparallel conventional vortices 
[HHIET], between coherence vortex-antivortex pairs. Non-equilibrium coherence simplex 
dynamics may reveal VCS reactions such as VCS intercommutation or reconnection 
events and coherence-simplex turbulence (cf. figure 4). 

In contrast to the topological arguments given earlier, a study of the dynamics 
of coherence simplices presupposes knowledge of the equations of motion underlying a 
given field ^/'(r, t). Suppose that this equation, for the field supporting a given coherence 
simplex or simplices, is given by the Heisenberg equation of motion 



^fi^^(r,t) = [V^(r,t),^], (9) 

where H is the Hamiltonian operator. From here on we focus on a free-field case 
(interactions alter the dynamics of coherence simplices but do not invalidate the 
topological considerations) H = —h'^V'^ /2m and define a linear differential operator 

in terms of which ^ may be expressed as 

£V^(r,t) = 0. (11) 

Acting with the components A; of £ on the coherence function g'^'P\xi,X2 ■ ■ ■X2p), we 
obtain a set of Wolf-like equations [TBI IIZl EH] 

Ckg^''\xi,X2...X2p) = 0. (12) 

These equations can be cast in a 'hydro dynamic' form [33] by separating the real and 
imaginary parts which yields two sets of coupled equations 

Reb(^)*£,(?(P)] = 0, (13) 
Im[^(^)*£fc(7(p)] = 0. 

These correspond to an Euler-like equation and a continuity equation for coherence 
flux, which govern the evolution of coherence in the system and hence the dynamics 
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of the coherence simphces which move with the coherence flow [25]. The spectrum 
of coherence excitations or coherence quasi-particles can be found by hnearising these 
equations [351 EZ] • 



The formal solution to (|12|) is 



l'^P\xi,X2 . . .X2p) = J J ■ J dx[dx2 . . . dx'^p (14) 

X G{xi, X2 ■ ■ ■ X2p, x^, X2 ■ ■ ■ X2p)g^'^\xi, X2 ■ ■ ■ X2p), 

which evolves the coherence function into the future using the forward-time many-body 
propagator 

G{Xi,X2 ■ ■ .X2p,x[,x'2 . . . X2p) = (15) 
- i{[lp{Xi) . . . ij{x2p)ij\x[) . . . ^Hx2p)])- 

Here we assume a time-ordering in which all primed time coordinates lie in the future 



of their corresponding unprimed coordinates. Equation (15) may be viewed as a 
generalized Huygens-type formula for the propagation of coherence [IH] , the propagation 
occurring from one set of equal-time space-time points {x[,X2, . . . ,X2p), to a different 
set of equal-time future space-time points (xi, X2, ■ ■ ■ , X2p)- For the special case of two- 
point correlation functions in one spatial dimension, this may be visualized as shown 
in figure 3(a). This depicts the propagation of the two-point equal-time one-spatial- 
dimensional correlation function g^^\ from its boundary value over the hyperplane 
ti = t2 = t, to its boundary value over the hyperplane ti = t2 = t + t, t > 0. Spatial 
coordinates are denoted with uppercase X, to distinguish from space-time coordinates 
denoted by a lowercase x. In this low-dimensional example, the coherence function 
over an infinitesimal patch ai A = {X[,ti = t,X2,t2 = t) = {x'i,X2) with infinitesimal 
area dX[dX2 = dx'^dx'2 in the X[ — X'2 plane (denoted Hi), propagates to give the 
contribution dx[dx2G{xi, X2,x'i, X2)g^^\x[, X2) to the coherence over the entire Xi — X2 
plane (denoted 112). This contribution is the coherence-theory analogue of the concept 
of a Huygens elementary wavelet, as embodied in the Huygens-Fresnel principle [16] . 
Every such contribution (from the points A,B,C etc. in figure 3(a)) is then summed, 
via Eq.(15), to give the propagated coherence function g^^\xi,X2) over the plane 112. 



This superposition of elementary Huygens-type wavelets may lead to points D at 
which g^^^ vanishes, due to what may be termed 'complete destructive interference of 
the coherence-function wavelets' at such a point. Coherence vortices (A/i and J\f2 in 
figure |4|^b), which appear as points F, G, L and M in the planes Hi and H2), for which a 
zero accompanied with a non-zero value of special case of such complete 

destructive interference of coherence. More 'regular' points such as E in figure 3(a) (and 
H, I, J and K in figure [4]^b)) correspond to non-zero values for g^^\ 

By encircling the coherence vortex G in figure |4|^b), the associated simple smooth 
closed path sampled at points H, I, J and K picks up a phase factor which identifies a 
vortical coherence field. If such a path integral yields non-zero ki, even if the field is 
fully or partially coherent at points H, I, J and K, there must be at least one point such 
as G in the space enclosed by the path, where coherence vanishes. A low-dimensional 
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Figure 3. Coherence simplices illustrated. As discussed in the main text, panel 
(a) shows a Huygens-type construction for the propagation of two-point correlation 
functions from hyperplane Hi to hyperplane 112. Three Huygens wavelets are shown, 
emanating from points A, B and C; r is a control parameter, which can be pictured 
as a generalized propagation distance, whose increasing value continuously evolves Hi 
into 112. Panel (b) shows an open nodal line Afi which threads coherence vortices 
located at points F and G, with A/2 denoting a closed coherence-vortex loop. A third 
possibility, that of a nodal-line knot A/a or a coherence Hopfion (cf. pT j [T2 t [T3 t fT4 l [T5] ) . 
is not shown. 



example, of this, is studied in Ref . [31] . Generically, the hypercore of a coherence simplex 
corresponds to a zero in the function g^^\xi, X2 ■ ■ ■ X2p) with an associated non-zero 
coherence circulation Kp. Note that the propagator, (16), being a correlation function 
itself, may also be vortical. 

Figure 3(b) also illustrates the conservation of the coherence circulation. Consider 
a coherence vortex A/i which at time t pierces the X[ — X'2 plane and is the only 
coherence vortex present at that time. As for conventional vortical systems [H [55] . 
the total circulation is a topological invariant and therefore must be preserved during 
the propagation of the coherence over any time interval r. Note, however, that 
coherence vortex-antivortex pairs (L, M) may nucleate and/or annihilate in the course 
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Figure 4. Vector coherence dynamics. The nodal Hnes corresponding to the 
propagation of the coherence field shown in figure 2a. Due to the symmetry of the 
problem only the (f, f) spin component of the coherence function is shown. Green 
and black lines correspond to m = 1 and m = — 1 coherence circulations, respectively. 
Coherence vortex loops and coherence vortex displacement waves are clearly visible in 
the figure. Coherence phase function at time i = 5 is also plotted (coloring corresponds 
to that in the figure 2 (b) and (d)). 



of propagation without affecting the value of the total coherence circulation. Such 
coherence vortex pairs appear as a manifestation of coherence vortex space-time loops, 
such as N'2, along which the coherence vanishes. Coherence- vortex nodal- line knots (not 
shown) are also possible; such knotted vortices have been studied both theoretically and 
experimentally, for the case of coherent fields [TT], |T21 HSl [IH US]- 

Finally, we demonstrate the evolution of vector coherence simplices, shown in figure 
2 (a,b) and confined in a toroid of length 12a, by propagating the coherence function 
forward in time as described by (15). The nodal lines of coherence with m = ±1 are 
denoted by green and black lines. Coherence vortex waves and loops of the kind A/2 (the 
fundamental elements of quantum coherence turbulence) are clearly visible in the figure. 
Space-time points where green and black lines join correspond to coherence vortex pair 
creation and annihilation events. 
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6. Discussion 

We have introduced a hierarchy of coherence simplices. A zeroth-order coherence 
simplex is nothing but a conventional quantized vortex in complex wavefunctions 
describing e.g. coherent optical fields or superfiuid matter. Higher-order vortical 
coherence simplices emerge as quantized vortices in generic multi-dimensional 
correlation functions. We have discussed the topological structure and dynamics 
of such coherence simplices. A hypercore of a coherence simplex corresponds to a 
topologically inevitable total destructive interference of the coherence-function wavelets. 
The existence of first-order vortical coherence simplices have already been verified 
experimentally using interferometric measurements of light |2l]. To our knowledge, 
they are yet to be discovered in the quasi-particle and matter-wave counterparts of 
coherent light. 

Furthermore, the Berezinskii-Kosterlitz-Thouless mechanism, which enables quasi- 
long-range order and superfiuidity in two-dimensional systems, relies upon the coherence 
of vortex-antivortex pairs [591 EO]- However, direct observation of such spontaneously 
forming vortex dipoles has remained elusive [HH E21 E3, EH ES]- It would be interesting 
to apply the theory of coherence simplices to such systems and to investigate the role 
of VCSs in the superfiuid-normal phase transition in low-dimensional systems. The 
recently created photon Bose-Einstein condensates [66] provide a natural platform to 
study coherence simplices and their relation to conventional quantized vortices which 
are a hallmark of superfluidity and are responsible for many of their special properties. 
It is therefore fascinating to contemplate the analogous idea that higher-order coherence 
simplices might underlie and signify some fundamental yet so far undiscovered physical 
properties of matter. Finally, it was pointed out in [31] that coherence simplices may 
be associated with the process of quantum mechanical decoherence — a suggestion which 
warrants further research. 
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